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Abstract 


We propose a class of actions for the spacetime metric that introduce cor¬ 
rections to the Einstein-Hilbert Lagrangian depending on the logarithm of some 
curvature scalars. We show that for some choices of these invariants the models 
are ghost free and modify Newtonian gravity below a characteristic acceleration 
scale given by ao = c/j, , where c is the speed of light and n is a parameter of 
the model that also determines the late-time Hubble constant: Hq ~ fi. In these 
models, besides the massless spin two graviton, there is a scalar excitation of the 
spacetime metric whose mass depends on the background curvature. This de¬ 
pendence is such that this scalar, although almost massless in vacuum, becomes 
massive and effectively decouples when one gets close to any source and we re¬ 
cover an acceptable weak field limit at short distances. There is also a (classical) 
“running” of Newton’s constant with the distance to the sources and gravity is 
easily enhanced at large distances by a large ratio. We comment on the possi¬ 
bility of building a model with a MOND-like Newtonian limit that could explain 
the rotation curves of galaxies without introducing Dark Matter using this kind 
of actions. We also explore briefly the characteristic gravitational phenomenol¬ 
ogy that these models imply: besides a long distance modification of gravity they 
also predict deviations from Newton’s law at short distances. This short distance 
scale depends on the local background curvature of spacetime, and we find that 
for experiments on the Earth surface it is of order ~ 0.1mm, while this distance 
would be bigger in space where the local curvature is significantly lower. 
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1 Introduction 


The relative importance of the gravitational interaction increases as we consider larger 
scales, and it is at the largest scales that we can measure where the observed gravita¬ 
tional phenomena do not agree with our expectations. The Hubble constant measuring 
the rate of expansion of the Universe does not fall with time as predicted by General 
Relativity (GR) for a Universe that contains only known forms of matter, and the 
dynamics of galaxies seem to require much more matter than observed if explained in 
terms of GR. The most common approach to these problems is to assume the presence 
of unseen forms of energy that bring into agreement the observed phenomena with 
GR. The standard scenario to explain the dynamics of galaxies consists in the intro¬ 
duction of an extra weakly interacting massive particle, the so-called Cold Dark Matter 
(CDM), that clusters at the scales of galaxies and provides the required gravitational 
pull to hold them together. The explanation of the observed expansion of the universe 
requires however the introduction of a more exotic form of energy, not associated with 
any form of matter but associated with the existence of space-time itself: vacuum en¬ 
ergy. And while CDM can be regarded as a natural possibility given our knowledge of 
elementary particle theory, the existence of a non-zero but very small vacuum energy 
remains an unsolved puzzle for our high-energy understanding of physics. However, the 
apparent naturalness of the CDM hypothesis finds also problems when one descends 
to the details of the observations. Increasingly precise simulations of galaxy formation 
and evolution, although relatively successful in broad terms, show well-known features 
that seem at odds with their real counterparts, the most prominent of which might be 
the “cuspy core” problem and the over-abundance of substructure seen in the simula¬ 
tions (see e.g. |T]). But, despite of this, the main problem that the CDM hypothesis 
faces is probably to explain the correlations of the relative abundances of dark and 
luminous matter that seem to hold in a very diverse set of astrophysical objects ,2J. 
These correlations are exemplified in the Tully-Fisher law j3j and can be interpreted 
as pointing to an underlying acceleration scale, below which the Newtonian potential 
changes and gravity becomes stronger. This is the basic idea of MOND (Modified 
Newtonian Dynamics), a very successful phenomenological modification of Newton’s 
potential proposed in 1983 4) whose predictions for the rotation curves of spiral galax¬ 
ies have been realised with increasing accuracy as the quality of the data has improved 
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|5]. Interestingly, the critical acceleration required by the data is of order ao ~ cH 0 
where Hq is today’s Hubble constant and c the speed of light (that we will set to 1 from 
now on). The problem with this idea is that MOND is just a modification of Newton’s 
potential so it remains silent in any situation in which relativistic effects are important. 
Efforts have been made to obtain MONDian phenomenology in a relativistic generally 
covariant theory by including other fields in the action with suitable couplings to the 
spacetime metric |6J (see also [7 for other approaches to galactic dynamics without 
Dark Matter). But these models do not address in a unified way the Dark Energy 
and Dark Matter problems, while a common origin is suggested by the observed co¬ 
incidence between the critical acceleration scale and the Dark Energy density. In this 
paper we will propose a class of generally covariant actions, built only with the metric, 
that have the right properties to address these problems in a unified way, and where 
the relation ao ~ Hq finds a natural explanation. The theories we will consider modify 
gravity in the infrared, making it stronger below a characteristic acceleration scale, but 
this is not their only characteristic feature. When we are in a situation in which the 
dominant gravitational field is external, like in table-top experiments on Earth (that 
measure the gravitational held of some probes embedded in the dominant background 
gravitational held of the Earth), we can also expect short distance modifications of 
Newton’s potential. 

Regarding the long distance modifications, the source-dependent characteristic dis¬ 
tance beyond which Newtonian gravity is modified in these theories, that we shall call 
r c , is given by 

GnM _ p 

9 o" ’ \-V 

r 2 2 

where G^ is Newton’s constant, M the mass of the source and /i is a parameter of 
the model that also determines the late-time Hubble constant: H 0 ~ p. This makes 
these models promising candidates to build a theory with a MOND-like Newtonian 
limit that could address the dynamics of galaxies without the need for Dark Matter. 
But when measuring the gravitational attraction between two probes in the external 
dominant gravitational field of a massive object of mass M, at a distance from its 
centre, we can also expect short distance modifications of Newton’s law for distances 
smaller than 
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If we plug in this expression the radius and mass of the Earth (with // ~ Ho), we get that 
for table top tests of Newton’s law performed on the Earth surface tsd rs./ 10 2 * cm. This 
range is very interesting because it is the range currently being probed by experiments 
jHj. But notice that the phenomenology of these theories is very different from the 
one expected from other theoretical considerations that also suggest a deviation from 
Newton’s law at that scale motivated by the cosmological constant problem 1 or the 
gauge hierarchy problem US]. The fact that this scale is the same in both cases is 
just a numerical coincidence. If we performed the same experiment on space, in the 
neighbourhood of the Earth’s orbit for instance where the dominant gravitational field 
is that of the Sun, the relevant mass and distance we should use in the previous 
estimation is the Sun’s mass and the Sun-Earth distance. In this case the “short 
distance” corrections would be expected in our theory at distances less than about 
10 4 m! This however does not mean that there should be big modifications to the 
motion of the planets or other celestial bodies. When the gravitational held we are 
measuring is that of the Sun, the corrections are suppressed at distances less than r c 
that is in this case of the order of lQpAU ~ 10 n A;m. 

These characteristic experimental signatures arise in our theory because of the 
presence of an extra scalar excitation of the spacetime metric besides the massless spin 
two graviton. But while the graviton remains massless, the extra scalar has a mass that 
depends on the background curvature. This dependence is such that for the models 
that we will be interested on, those that modify gravity at large distances, this field 
becomes massive and effectively decouples when the background curvature is large, and 
in particular when we approach any source. 

In the next section we will briefly review the results we obtained in mm studying 
models that involve inverse powers of the curvature in the action, giving some gen¬ 
eral expressions and discussing the generic features of the framework we will use for 
modifying gravity in the infrared. In particular we will focus on a class of models 
that had been proposed to address the acceleration of the Universe ca and also mod¬ 
ify gravity at large distances mm • This discussion will enable us to motivate the 

1 It is well known that naturalness arguments lead to the expectation that new physics associated 

with electro-weak symmetry breaking, besides the Higgs boson, should be seen in the LHC. Otherwise 

the electroweak scale becomes unstable under quantum corrections. Applying the same logic to the 

gravitational sector one would expect new gravitational phenomena to kick in at the vacuum energy 
scale that would cut-off the quantum divergences contributing to the vacuum energy. This hypothetical 
new physics should be seen in sub-mm measurements of Newton’s potential (5j ■ 
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class of actions that we will present in the third section, that depend on the logarithm 
of some curvature invariants. We will see that the theories that we propose in this 
section modify gravity at the MOND characteristic acceleration scale, and the gravita¬ 
tional interaction can easily become stronger at large distances. We will explore briefly 
the characteristic gravitational phenomenology expected in these models, discussing 
possible tests of these theories. In the fourth section we offer the conclusions. We 
will comment on further generalisations of the proposed actions and on the generic 
phenomenological features expected in the class of models that modify gravity at the 
MOND acceleration scale. We will also comment on the possibility of obtaining these 
theories as an effective action for the spacetime metric that takes into account strong 
renormalisation effects in the infrared that might appear in GR. 

2 Modified gravity as an alternative to Dark En¬ 
ergy 

Recently, models involving inverse powers of the curvature have been proposed as an 
alternative to Dark Energy nam- In these models one generically has more propagat¬ 
ing degrees of freedom in the gravitational sector than the two contained in the massless 
graviton in GR. The simplest models of this kind add inverse powers of the scalar cur¬ 
vature to the action (A£ cx 1/R n ), thereby introducing a new scalar excitation in the 
spectrum. For the values of the parameters required to explain the acceleration of the 
Universe this scalar field is almost massless in vacuum and one might worry about the 
presence of a new force contradicting Solar System experiments. However in a recent 
publication cu we showed that models that involve inverse powers of other invariants, 
in particular those that diverge for r —> 0 in the Schwarzschild solution, generically 
recover an acceptable weak field limit at short distances from sources by means of a 
screening or shielding of the extra degrees of freedom at short distances 

But let us start by discussing the linearisation and vacuum excitations obtained 
from generic actions built with the Ricci scalar and the scalars 

P = R IW R IW and Q = R^ Xp R^ xp . (3) 

If the Lagrangian is a generic function C = F(R, P , Q) the equations of motion for 
the metric will be of fourth order and we can expect that the particle content of the 
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theory will have eight degrees of freedom: two for the massless graviton, one in a scalar 
excitation and five in a ghost-like massive spin two field ra- Expanding the action 
in powers of the curvature perturbations it can be seen that at the bilinear level the 
linearisation of the theory over a maximally symmetric spacetime will be the same as 
that obtained from fT TT^ 


S = Jd^xy^g-^ 


-A + 5R + 


6 ml 


R 2 


I67 tGn 

where C^ u \ a is the Weyl tensor and we have defined 
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A = (F — RF r + R 2 (F rr /2 — Fp/A — Fq/6) + R 3 (F RP /2 + F R q/3) 


+R* (Fpp/8 + Fqq/18 + F pq /6)) 0 ( 5 ) 

<5 = (Fr — RF rr — R 2 (F rp + 2F R q/3) 

-R 3 {Fpp/4: + Fqq/9 + F P q/3)) 0 ( 6 ) 

m 0 2 = (( 3 F rr + 2F P + 2Fq) + R ( 3 F RP + 2 F R q) 

+R 2 (3Fp P /4 + Fqq/3 + F pq )) 0 (7) 

mF = — (F P + 4Fq) 0 . (8) 


Here < ... >o denotes the value of the corresponding quantity on the background and 
Fji = d R F } etc... One can see that the situation for the perturbations over vacuum 
in any modified theory of gravity (built with R , P and Q ) will be the same as in 
Einstein gravity supplemented with curvature squared terms. It is well known that for 
the action (ED the mass of the ghost is ~ 7n 2 and that of the scalar is ~ m$. So in 
the case in which F(R, P, Q ) = F(R, Q — 4 P), m^ 2 = 0 and there is no ghost in the 
spectrum, but there is still the extra scalar. It is easy to check that in the models that 
explain the acceleration of the Universe by using actions that involve inverse powers 
of the scalar curvature alone, the mass of the scalar is proportional to some positive 
power of the scalar curvature era And notice that including other terms in F and 
fine-tuning the parameters we can make 2 = 0 in vacuum, as suggested in EE], by 
including R 2 corrections to the action besides the 1/Rot Log(f?) ones. But the infinite 
mass (or absence) of the scalar is a property only of a particular background in these 
models. If we evaluate the expression for m o for F = R — /i 4 /R + R 2 /m 2 for instance, 
we get 

/ ,.4 1 \ 
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and choosing a particular value of m s , namely m 2 = 3\/3fi 2 , we can make mo —»► 
in vacuum, where (R) 0 = \/3r 2 . But we see that when the scalar curvature is bigger 
than its vacuum value, the scalar mass returns to its natural value, mo ~ m s ~ //, 
which shows that in these models the scalar is still present in general, and its mass is 
very small in most situations 2 . But for the models that include only inverse powers of 
the curvature, besides the Einstein-Hilbert term, it is however possible that in regions 
where the curvature is large the scalar has naturally a large mass and this could make 
the dynamics to be similar to those of GR ra. But the scalar curvature, although 
bigger than its mean cosmological value, is still very small in the Solar System for 
instance. So, although a rigorous quantitative analysis of the predictions of these 
models for observations at the Solar System level is still lacking in the literature, it is 
not clear that these models constitute a viable alternative to Dark Energy because one 
can expect that the effects of this extra held should have been observed 3 . 

But the story is different if we include inverse powers of curvature invariants, like 
Q , that grow at short distances in the Schwarzschild solution. In this case the mass of 
the extra scalar held is guaranteed to grow as we approach any source and we recover 
Einstein gravity at short distances uma- For studying these effects we focused on 
actions of the type 

r i r u An+2 

S = Jd x V^9 WnGN R ~ ( ai?2 _ 4 P + Q)n 

These actions, with /i of the order of the late time Hubble constant Hq, were proposed 
hi EE] motivated by the Dark Energy problem and have cosmological solutions pro¬ 
viding a good ht to the SN data [2T]. When linearising over vacuum we hnd, besides 
the usual massless graviton, a scalar excitation of the metric with a mass of order 
mo ~ n ~ Hq. But in the expanded action higher order non-renormalisable operators 
appear suppressed by inverse powers of the background curvature. This means that the 
linearisation will break down when the energy of the fluctuations is extremely small: 
the strong coupling scale in this theory is A s ~ (jSd p Hl n+A /ji 2n+r ) ~ ( M p H q) 1 ^ 4 T?j . 
This, in turn, means that for a spherically symmetric solution, the linearisation over 

2 This is in contrast with the absence of the ghost for actions in which P(P, P, Q) = F(R, Q — 4P). 
In this case m^ 2 = 0 independently of the background curvature. 

3 There is an alternative formulation of these models in which the connection and the metric are 
varied independently, the so-called Palatini formalism, and in this case these theories might be viable 

EDI- 
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vacuum will break down at a huge distance from any source, at the so-called Vainshtein 
radius, that is for these theories ry ~ (GnM /i 2n+1 /// ( 2n+4 ) ~ [GnM/ Hq) 1 ^ 4 |T2] . At 
smaller distances we can not trust the results obtained using this expansion. It is also 
important to keep in mind that in this long-distance regime, where we can apply the 
linearisation of the action over vacuum, the effective Planck mass is given by (see t I2|) 


M 2 , ,,, = M 2 

P(eff) P 


-1 


<5 + 


Wo 

3?tIq 


= 1 + 


6 n(a — 1) 

[n + l)(6a — 5) 


Ml 


( 11 ) 


where M 2 = (87tGat) 

As we have said, at distances less than 7y we enter a non-perturbative regime 
where we can no longer use the linearised action over vacuum because higher order 
non-renormalisable operators become more important than those involving only two 
powers of the fluctuations. But we can get some information about the short dis¬ 
tance behaviour of the solutions by noticing that the extra term that the modification 
introduces will be unimportant at short distances in the spacetime of a spherically 
symmetric mass. The reason is that this term will always be suppressed by inverse 
powers of Q, that for the Schwarzschild solution reads Q = 48(G nM) 2 /r 6 and grows 
at short distances. This tells us that any spherically symmetric solution for which the 
curvature grows at short radius will converge to the Schwarzschild one with a Planck 
mass given simply by M 2 = (87rGA) _1 - We can then make a different kind of expan¬ 
sion, a weak held expansion over the Schwarzschild solution that shows that in these 
models Newtonian gravity is modified beyond a distance given by EU 


3n+2 = 

1 c 


(G N M) n+l 

H 2n+1 


( 12 ) 


In fact, for these theories, the spherically symmetric solution, in an expansion in powers 
of r/r c reads: 


ds 2 




+ 


2G n M 


a 


6n+4 


+ 0((r/r c ) 12n+8 ) 


dt 2 


(13) 


2GnM ( 1 + a ^ 6n 2 + 3 ) (l) “ + o ((r/r^ 12 ^ 


„) 12n+8 ) 


dr 2 + r 2 dVt\ , 


where M is the mass of the object sourcing the held and a = ( 6 r l+ 3 ylll 3 n. ■ 

The distance r c is very large, at least of the order of parsecs for a star like the 
Sun (assuming n > 1), so from this expansion we see that the corrections that we 
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Figure 1: Schematic representation of the structure of the theory on the en¬ 
ergy/distance scale. 


can expect in these theories at the Solar System level are very small. Unfortunately, 
this short distance expansion breaks down for distances of order r c , and notice that 
r c -C ry. So in these theories we are led to the following picture: at small energies 
(E < A s ) or large distances (r > ry) we have a scalar tensor theory, with an almost 
massless scalar and an effective Planck mass given by (ED- As we increase the energy 
of the fluctuations or go to shorter distances we enter a non-perturbative phase, but at 
even higher energies (E > Aqr) or smaller distances (r < r c ) we can neglect the effects 
of the modification and we recover the standard GR dynamics, with a Planck mass 
given by A4 P . Here A gr = (M” +1 /r 2n+1 ) 1 A 3n+ A the energy scale associated with this 
high energy recovery of Einstein gravity. This situation is depicted in hg.l. Since all 
the expansions that we have used break down in this non-perturbative intermediate 
regime, at this point we can only assume that the dynamics are consistent for this 
range of energies/distances, and that one can obtain a consistent matching between 
the long distance and short distance solutions. However we regard this as a model 
dependent issue, and we can extract already a lot of useful information just knowing 
the dynamics in the low and high energy regimes. 

In terms of the particle content of the linearised theory we can get some “non- 
perturbative” insight in the reasons behind the behaviour of these solutions by evalu¬ 
ating the expression for mo in them. We see then that we can expect that the mass of 
the extra degree of freedom that the modification introduces has a contribution in the 







spacetime of a spherically symmetric mass that goes like t I2j 

, ,, Qz_ (GnMY^_ 

& Source™ sif) ~ ^2n+l ~ ^2n+l r 3n+3' ( ' 

This effective mass decouples the scalar at short distances, where Q /i 4 , and we 
recover Einstein gravity in this domain. It is also worth to point out that one recovers 
the distance r c as the distance at which m s {r) ~ r _1 , with m s (r) given by the expression 
above. Since m s {r ) grows faster than r _1 , at smaller distances the scalar effectively 
decouples. But the r-dependent mass of the scalar field is not the only effect of the 
modification: as we have seen at large distances from sources the effective Planck mass 
is given by eq. m while at short distances it is simply M p . So there is a rescaling or 
“running” of the Planck mass with the distance to the sources. 

However, one can see that these theories, although capable of explaining the acceler¬ 
ation of the Universe without Dark Energy, are not capable of explaining the dynamics 
of galaxies without Dark Matter. First, the critical distance at which the modification 
becomes important (r c ) is too large, since the data consistently indicate that any such 
theory should have noticeable effects at a distance r c given by eq.(P). And second, 
the effective Planck mass that we obtain at large distances is not significantly reduced 
with respect to the one that we get at short distances except for a very small range of 
values of a that, at least for the n = 1 case, are not phenomenologically acceptable EH 
We need the reduction of the effective Planck mass in vacuum in order to get a large 
enhancement of the gravitational interaction at large distances as required to have any 
hope to fit the data without introducing Dark Matter. The critical distance that we 
obtain in these theories does however suggest that in the n —> 0 limit the modification 
becomes important at a distance that corresponds to the MOND characteristic accel¬ 
eration. This motivates the use of Logarithmic actions as a possibility for obtaining 
a modification of gravity as an alternative to Dark Matter. A first assessment of this 
possibility is the goal of the next section. 

3 Logarithmic actions and modified gravity as an 
alternative to Dark Matter 

As we said, the considerations of the previous section motivate the introduction of an 
action depending on the logarithm of the curvature in order to get a modification of 
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gravity at the MOND characteristic acceleration scale. In this section we will discuss 
what appears to be the simplest possibility, namely actions of the type 

S = J {R- Ai 2 Log [f(R, Q-AP)]} , (15) 

where for the function / we will only assume that 

/-0 for R'jtvx —> 0 , (16) 

and we can approximate 

/~Q/Qo when Q^>R 2 ,P. (17) 

In this case Minkowski spacetime will not be a solution of the theory but there will 
typically exist de Sitter solutions with H 0 ~ /i. Notice that the addition of a cosmo¬ 
logical constant will not change the form of the action 4 . As we will see we can apply 
the discussion of the previous section to these theories simply taking n — 0 in the rel¬ 
evant formulae. We will study in the next subsection the behaviour of the spherically 
symmetric solutions of this theory at short distances from sources, and we will show 
that in this domain the corrections to the Schwarzschild geometry are small. Also, it 
will be made clear that the Newtonian potential has large corrections at a distance 
given by r c in eq. (HD- In subsection 3.2 we will study the linearisation of this theory 
in vacuum. We will give the conditions for the stability of de Sitter space and we will 
see that the effective Planck mass in vacuum can easily be reduced with respect to the 
one at short distances by a large ratio, and in this case gravity becomes significantly 
stronger at large distances. The expansion of the action in vacuum is the relevant one 
to apply in some dynamical situations at very large distances from sources and in late¬ 
time cosmology, but one should keep in mind that this linearisation breaks down when 
one approaches any source or when the ambient curvature is large (like, for instance, 
in the early Universe). In subsection 3.3 we will discuss some of the characteristic 
experimental signatures that these theories imply. 

3.1 Short distance solution 

In this subsection we will follow the same strategy as in m in order to study the 
behaviour of the solutions at short distances from sources. The spherically symmetric 

4 Adding a cosmological constant would simply redefine the function /. When writing the action 
in this form we are of course assuming that the mass scales appearing in / are “reasonable” in terms 
of n so that H 0 ~ /i. 
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solutions for this theory are obtained by solving the equations 


Gfiu + nu — 0 , 


( 18 ) 


where G^ is the usual Einstein tensor and /i 2 //,,,, is the extra term generated by the 
logarithmic part of the action. One can see that this extra term, when evaluated in 
the Schwarzschild solution, is subleading with respect to the terms that appear in the 
Einstein tensor when r <C r c . This indicates that as r —> 0 the corrections with respect 
to the Schwarzschild geometry will be small. We can then consider a small perturbation 
of the black hole geometry and solve at first order in the perturbations. So we take 
the ansatze 


ds 2 = — 


1 - 


2GnM 


+ eA(r) 


dt 2 + 


2 G n M 
1-b eB(r) 


i -i 


dr 2 + r^dn '2 , (19) 


and treat e as a small expansion parameter. We can expand the full equations in powers 
of e: 


G^u - G ( °J + eG^J + e 2 Gj?J + ... , 

H „„ = H<$ + eHj£+e*H<$ + .... ( 20 ) 

Since the Schwarzschild solution solves the ordinary Einstein equations, we have (j,',]! — 
0. Treating n 2 as an order e parameter, at first order in our expansion the equations 
for A and B become (from now on we set e = 1) 


g;'. ( 1) = -y k (0) ■ 


( 21 ) 


For the tt component of this equation we End: 


B + rB' 




r 2 2 G n M 

while the rr component reads: 

{2G N MA- rB){2G N M- r)- 1 + rA' p 2 


13GatM — 8 r + G^M Log 


48 {G N Mf 

GQ 0 




2G n M 


7G n M — 2 r — G n A/Log 


( 22 ) 


48(GivM) : 

r 6 Qo 

(23) 


We can solve the previous equations yielding 

4 


B(r) = 
A(r) = — 


2 G n M\ r 


r 


rv 


G n M 
3 r 


15 + Log 


2G n M\ (r 


4 r 


4 G n M 


3 r 


5 — Log 


48 {G N Mf 
Q 0 r e 

'48 {G N Mf 
Qor 6 


, (24) 
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where r c is now given by eq. (HD- So we see that at short distances the corrections to 
Newton’s potential are indeed suppressed by powers of (r/r c ) 4 . In fact, we can con¬ 
sider our approximate solution as the first order in the Taylor expansion in powers of 
r/r c of the exact solution, where the next order in the expansion would be O ((r/r c ) 8 ). 
At distances of order r c the expansion breaks down, and it is clear that we can ex¬ 
pect a significant modification of Newton’s potential for larger distances (or smaller 
accelerations). A difference with respect to the actions involving inverse powers of 
the curvature is that now the scalar curvature does not go to zero at short distances. 
This is so because the logarithmic part of the action does not go to zero for r —> 0, 
as happened in the case of the inverse powers in the previous section, but the scalar 
curvature (or the extra term in the action) is of the order of the one in vacuum so the 
effect is very small. 

In the next subsection, we study the linearisation of this theory in vacuum, which 
will enable us to obtain the ultra large distance behaviour of the spherically symmetric 
solution corresponding to a mass source. 


3.2 Linearisation in vacuum 


In vacuum, without any matter source, the action has de Sitter solutions with curvature 
R = 12Hq where Hq can be found as a solution of 

6 Hi = y ^Log/ - 6 Hi /b - 2 ° g ° /o \ . (25) 

To check the stability of these solutions one can consider the modified Friedmann 
equation in vacuum obtained in this theory: H = h(H, H) and expand for small 
perturbations H and H = H — H 0 of the equation of motion. We find: 

H « h H (H 0 , 0 )H + 0 )H = 16 Hi Q + H - 3 H 0 H , (26) 


where we have defined 


4 O' 2 

Ci = 5 H-and C 2 = 


m 


.is m 

ml 


(27) 


m 0 and S have been defined in the second section. One can easily check that the 
fixed point (H = 0, H = 0) of this dynamical system is an attractor (repcller) if the 
coefficient in front of H is negative (positive), so de Sitter space will be stable as long 


12 





as 

4 Ci , . 

1 + 7T 1 < 0. (28) 

One would obtain the same result by analysing the equation of motion of the prop¬ 
agating modes, and it is instructive to do so. Because the action is a function of Q 
and P only through the combination Q — 4 P, the ghost is absent and there is just 
a massive scalar held in the spectrum besides the massless graviton. So this held is 
the only possible source of instability. It was shown in jH2j that in this case de Sitter 
spacetime is stable as long as m 2 > — 9H^/A, and the mass of the scalar is 

< - -«0 2 (f + l«g) . (») 

which is consistent with the phase space analysis presented here and also with the 
results of ra, for F(R) theories. But even when de Sitter space is unstable one should 
not disregard the model. It has been shown that for the actions involving inverse powers 
of the curvature de Sitter space is unstable in many cases but the late time background 
corresponds to a power law FRW cosmology that is nevertheless phenomenologically 

interesting D3E5- 

An important feature of this linearisation is that the effective Planck mass that we 
obtain in vacuum controlling the coupling of the spin two graviton, is now given by 

=Ci M i =(i - y - ~ 2 f R!v ) K (3°) 

At short distances however the value of the Planck mass is just M p , as we saw using 
the short distance expansion over Schwarzschild geometry. But in vacuum R 2 ~ P ~ 
Q ~ yU 4 ~ Hq and we can expect that C\ will depart significantly from 1. This means 
that when this linearisation is applicable, we get a rescaled Planck mass with respect 
to the one we would infer at short distances and one can expect an enhancement or 
suppression of the gravitational interaction at large distances. In fact, applying this 
linearisation, we get for a spherically symmetric mass the solution 

r 2 ^ dr 2 + r 2 dfll (31) 


ds 2 ~ 


8 G% ff) M 
3 r 


H, 


2„2 


dt 2 + 1 


3 r 


Hi 


where G— Gn/C\ is the effective long-distance Newton’s constant. Remember 
that the linearisation we are using to get this result breaks down at a huge distance 
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from any source, at the Vainshtein radius ry ~ (GArM/i/g) 1 / 4 . At shorter distances 
we can not use the linearised version of the theory, and the solution above is only a 
good approximation for r > ry. But we have seen that at even shorter distances, 
r < r c , we can use a different expansion: the one that we considered in the previous 
sub-section that shows that in this regime the solution, eas. (|19l24j) . approaches the 
Schwarzschild one with a Planck mass given by . So also in these theories there is 
a non-perturbative regime in the range of distances r c < r < ry, where one should get 
an interpolation between the solutions eqs. (USED and eq. (ED). In this intermediate 
range of distances both the linearisation of the theory over the vacuum and the short 
distance expansion over the Schwarzschild solution break down, but we will assume 
that the function f(R , Q — AP) is such that a consistent matching exists. When C\ < 1 
and Newton’s constant is enhanced at large distances, the interpolating regime would 
play the role of a “dark halo” surrounding any source. 

Even though we have not solved the gravitational equations in this non-perturbative 
regime, we would like to stress some important characteristics of this intermediate 
range of distances that make these theories promising candidates to build a MOND-like 
modification of gravity that could explain the dynamics of galaxies without recourse 
to Dark Matter. First, the modification becomes important below a characteristic 
Newtonian acceleration scale determined by a parameter (/r) that also determines the 
late-time Hubble constant, as consistently indicated by the data [5jJ. And second, the 
gravitational interaction can be enhanced at long distances by a ratio that depends 
on the particular function / that we choose, but that can easily be of the required 
magnitude. As an example we plot in fig.2 the ratio = 1/Ci, for / = 

(2 R 2 — 4 P + Q)/Q o, as a function of Log[/r 4 /Qo]/2. The apparent mass discrepancy 
that one would infer cosmologically or at large distances is of the order of this ratio. 
There are further corrections because the theory at large distances is of the scalar-tensor 
type, with an almost massless scalar, so we get an extra 4/3 factor in the Newtonian 
potential and the deflection of light would also get an analogous correction factor. 

These theories make a diverse variety of characteristic experimental predictions that 
can be used to test or falsify them. In this section we have seen the approximate form 
of the solution for a spherically symmetric mass at short distances (r < r c ) and at 
long distances (r > r v ) in de Sitter space. Using these approximate solutions and the 
linearised version of the theory, once we have chosen a particular function /, we can al- 
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Figure 2: Ratio of the effective Planck mass at short distances (r < r c ) over the one at 
large distances (r > ry) for a spherical mass in the model (11511 with / = (2i? 2 — 4P + 
Q)/Qo, as a function of Log /Q l J 2 }. This represents roughly the enhancement of the 
gravitational interaction at large distances. 


ready compare the predictions of the theory with precision Solar System measurements 
and a wealth of cosmological, astrophysical and large scale structure data. But these 
are not the only testing grounds for the theory. It is also possible to obtain predictions 
that would differ from those of GR for laboratory experiments measuring the gravita¬ 
tional interaction of small bodies. The discussion of these experimental implications is 
the subject of the following subsection. 

3.3 Experimental implications 

The theories we are studying in this section offer a diverse range of characteristic 
experimental predictions differing from those of GR that would allow their falsification. 
The most obvious tests would come from the comparison of the predictions of the theory 
to astrophysical and cosmological observations where the dynamics are dominated by 
very small gravitational fields. But we would have also some small effects for the 
motion of the planets or other celestial bodies in the Solar System and as we have said 
short distance modifications of Newton’s law. In the following we will briefly discuss 
separately these possible tests, and we will give order-of-magnitude estimations of the 
expected effects. 
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3.3.1 Short distance deviations from Newton’s law 


Probably the most characteristic experimental signature of our theory corresponds 
to deviations from Newton’s law that should be seen at short distances, where the 
precise distance at which the modification is noticeable depends on the local background 
curvature. While a detailed computation of these corrections will be deferred to a future 
publication |2Bj, to see why this is so we can use the following arguments. As we have 
said, for theories of the type (USD there is an extra scalar excitation of the spacetime 
metric besides the massless graviton. The mass of this held is given in de Sitter space 
by eq.fEHD, but it picks up an r-dependent contribution in the spacetime of a spherically 
symmetric mass. In a generic background, if we study a system in a scale that reduces 
the space-time region of interest to be small enough, we can consider an effective theory 
in which the extra degree of freedom has a mass given by its local value. Now, if we 
evaluate the expression for mo in the spacetime of a spherically symmetric mass we get 
that at a distance of a massive object of mass M, the mass of the scalar is locally 
of order 


2 Q ( G N Mf 

rri — nu - 

s /i 2 r> 2 


(32) 


where we are assuming that Q i? 2 , P so / ~ Q/Qo- So in the effective gravitational 
theory that we should apply on the Earth surface there is, besides the massless spin two 
graviton, an extra scalar held with gravitational couplings and a mass given roughly by 
the expression above. A peculiar feature of this local effective theory on a Schwarzschild 
background is that there will be a preferred direction and this will be reflected in an 
anisotropy of the force that this scalar excitation will mediate [221 • But for the purposes 
of this section, an estimation of the expected order-of-magnitude of the corrections, we 
will simply focus on the value of the mass of the scalar. From effective held theory 
arguments we can then expect short distance modifications of Newton’s law, suppressed 
by tsd ~ miT 1 ~ (9(0.1 mm), when measuring the gravitational held of some probes 
on the Earth surface. This order of magnitude is very interesting because it is the one 
currently being probed by experiments [Hj. As we said there are also other motivations 
for expecting short distance modifications of gravity at this scale eng, but as we also 
said the coincidence of these scales is just a coincidence. If one was to measure the 
gravitational held of a small object in space, in the neighbourhood of the Earth’s orbit 
but far from the Earth, the relevant mass and distance we should apply in eq. fl22D is 
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the mass of the Sun and the Sun-Earth distance, since the Sun provides the dominant 
gravitational field in that situation. In this case we get mj 1 ~ 10 4 m. But notice 
that to observe a significant modification in the gravitational field of an object we 
have to measure it at a distance bigger than r c for that object, otherwise the self¬ 
shielding of the extra scalar excitation induced by the object itself is enough to switch 
off the modification. This means that locally in the inner Solar System we could only 
see significant modifications in the gravitational field of objects whose characteristic 
distance r c is smaller than 10 4 m, so its mass would have to be below ~ 10 9 kg. As 
an example we can mention that for an object of mass 10 3 kg orbiting the Sun at 
the same distance as the Earth, one can expect modifications of its gravitational field 
starting at a distance r c ~ 10 m (at shorter distances the scalar effectively decouples 
because of the gravitational field of the object itself) and extending to rsD ~ 10 4 m 
(at longer distances the mass induced by Sun’s gravitational field effectively decouples 
the scalar). 


3.3.2 Solar System observations 


Although the corrections that our theories introduce with respect to GR for a spheri¬ 
cally symmetric solution are suppressed by powers of (r/r c ) 4 , and this is a very small 
number for the Sun within the Solar System, the precision with which the motion of 
the celestial bodies of the Solar System is known makes it conceivable that one could 
see the corrections induced by this modification (similarly as happens in other long¬ 
distance modifications of gravity as the DGP model |24|h For instance, in the case of 
the precession of the perihelion of the planets, the anomalous shift (A</>) induced by a 
small correction (5V) to Newton’s potential (V)v) is given in radians per revolution by 
(see e.g. [21]) 


. , d ( 9 d 
A0 ~ 7Tr— r“ — 
dr \ dr 


SV 


so the correction induced by our modification increases with distance as 

4 

A(j) ~ 167T I — 

V r r 


(33) 


(34) 


From this expression we see that for the inner planets the perihelion shift induced by 
our modification is very small, it only becomes of the same order as the correction 
introduced by GR (with respect to Newtonian gravity) for Jupiter. But for the outer 
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planets the contribution to the perihelion shift of the modification can be dominant over 
the standard GR one 5 . Unfortunately the only reliable data regarding the planetary 
perihelion comes from the inner planets of the Solar System |2S], where our correction is 
negligible. But the best measured orbit is that of the Moon around the Earth by virtue 
of the Lunar Laser Ranging (see e.g. J2EI)- Using this the Earth-Moon distance is known 
with a precision of centimetres. Applying the formula above, for the Moon our theory 
predicts an anomalous shift of A <f> ~ lO -12 , to be compared with the achieved accuracy 
of 2.4 x 10 11 [21 • Although these numbers are not far from each other remember that 
we are just doing an order of magnitude estimation, since the parameter g is related 
to H 0 only after choosing a particular function /. So although these theories suggest 
the possibility of surprises in high precision astrometrical measurements, the numbers 
that we obtained also expose the difficulties of ruling out these theories using these 
effects, since for this we would need an improvement on the bounds of the anomalous 
precession of the Moon of at least two orders of magnitude. 

3.3.3 Cosmology and Astrophysics 

In this subsection we will comment on the possible comparisons that one could make 
of the predictions of these theories to astrophysical and cosmological observations (see 
for instance El), although any actual fit to these data is beyond the scope of the 
current paper. In this respect, once a particular function / is chosen, one can make 
unambiguous predictions for the rotation curves of spiral galaxies with the mass-to- 
light ratio being the only free parameter, since the normalisation of g will be fixed by 
the Hubble constant, and in fact this issue has been our main motivation for proposing 
actions of the type (USD- Our result for Newton’s potential at short distances (El can 
be parameterised as: 

V(r) = —v(x) , (35) 

r 

where x = r/r c and v(x) ~ 1 + |a; 4 + ... for small x. The challenge now is to find a 
form for / that yields a MOND-like phenomenology in the intermediate region. More 
specifically we need that v(x) ~ —x ln(x) for large x. The potential would then give rise 
to flat rotation curves that obey the Tully-Fisher law |3j. But also other aspects of the 
observations of galactic dynamics can be used to constrain a MOND-like modification 

5 Not only our correction increases with distance, but the GR one (A< j> ~ QitGnM/ r) decreases for 
larger radius. 
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of Newton’s potential, see e.g. |28| . And notice also that our theory violates the 
strong equivalence principle, as expected for any relativistic theory for MOND Tj, 
since locally physics will intrinsically depend on the background gravitational field. If 
we consider for instance a local system like an open cluster, rotating in the background 
field of the Milky Way, the solution (ESI) that we obtained for an isolated system is not 
necessarily valid anymore. This will be the case if the background curvature dominates 
the curvature induced by the local system, similarly to the “external field effect” in 
MOND. Notice that in order to establish the relevance of the external field effect in 
this framework we have to compare the induced external and internal curvatures, and 
not the accelerations. This is relevant for the case of the globular clusters orbiting 
the Milky Way, where MOND (or DM) evidence has been found where the internal 
acceleration is below ao, even if the external acceleration due to the Milky Way is 
bigger than ao [2HJ- However, if one compares the values of the curvature (i.e., the 
invariant Q) one finds that the curvature is dominated by the internal one, justifying 
the neglect of the external field effect and the treatment of these globular clusters as 
isolated objects, into the MOND regime. So the globular cluster data of j2H] favors 
models, like the one presented in this paper, in which the relevance of the external field 
effect is determined comparing curvatures, and not accelerations. 

At larger scales, where one can use the equivalence with a scalar-tensor theory 
more reliably, we can expect from continuity arguments that such a form for / (if it 
exists) will enhance Newton’s constant. One can then compare the theory against the 
observations of gravitational lensing in clusters, the growth of large scale structure and 
the fluctuations of the CMB. In fact, it has been pointed out that if GR was modified 
at large distances, an inconsistency between the allowed regions of parameter space 
would show up for (non-modihed) Dark Energy models when comparing the bounds 
on these parameters obtained from CMB and large scale structure EH- This means that 
although some cosmological observables, like the expansion history of the Universe, can 
be indistinguishable in modified gravity and Dark Energy models, this degeneracy is 
broken when considering other cosmological observations and in particular the growth 
of large scale structure and the Integrated Sachs-Wolfe effect (ISW) have been shown 
to be good discriminators for models in which GR is modified [39. 

Regarding the ISW effect, we would like to mention some characteristics of these 
theories that point to the possibility of getting a suppression of the low multipoles of 
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the CMB with respect to a ACDM cosmology, a feature shown by the data that can 
not be easily accounted for in the ACDM model. It has been recently pointed out that 
the fact that in the DGP model the effective Newton’s constant increases at late times 
as the background curvature diminishes, causes a suppression of the ISW that brings 
the theory into better agreement with the CMB data than the ACDM model [32] - In 
our case we can expect an analogous effect, so that the effective Newton’s constant for 
the cosmic perturbations depends on the background curvature in such a way that it 
increases at late times as the curvature diminishes. But despite this potentially good 
features it remains to be seen if one can get a good fit to the CMB data in the absence 
of Dark Matter in these models. 

4 Conclusions 

In this paper, motivated by the phenomenological success of MOND fitting the rotation 
curves of spiral galaxies without requiring Dark Matter, we have proposed a class of 
actions that modify gravity below the characteristic acceleration scale required by 
MOND: do ~ Hq. There are two effects in these theories that are responsible for 
the infrared modification. First, there is an extra scalar excitation of the spacetime 
metric besides the massless graviton. The mass of this scalar field is of the order of the 
Hubble scale in vacuum, but its mass depends crucially on the background over which 
it propagates. This dependence is such that this excitation becomes more massive as 
we approach any source, and the extra degree of freedom decouples at short distances 
in the spacetime of a spherically symmetric mass. This feature makes this excitation 
to behave in a way that reminds of the chameleon field of j2|, but in our case this 
“chameleon” field is just a component of the spacetime metric coupled to the curvature. 
But there is a second effect in these theories: the Planck mass that controls the coupling 
strength of the massless graviton also undergoes a rescaling or “running” with the 
distance to the sources (or the background curvature). This phenomenon, although 
a purely classical one in our theory, is reminiscent of the quantum renormalisation 
group running of couplings. So one might wonder if actions of the type csd could be 
an effective classical description of strong renormalisation effects in the infrared that 
might appear in GR (see e.g. [331 anc I references therein), as happens in QCD. In fact, 
corrections depending on the logarithm of the renormalisation scale are ubiquitous in 
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quantum field theory, and it appears natural to identify the renormalisation scale with 
a function of the curvature if we want to build an effective classical action for the 
spacetime metric that takes into account these quantum effects. Indeed, we have seen 
that these models offer a phenomenology that seems well suited to describe an infrared 
strongly coupled phase of gravity: at high energies/curvatures we can use the GR action 
or its linearisation as a good approximation, but when going to low energies/curvatures 
we find a non-perturbative regime. At even lower energies/curvatures perturbation 
theory is again applicable, but the relevant theory is of scalar-tensor type in a de Sitter 
space. 

We would like also to emphasise that there are clearly many modifications of the 
proposed class of actions that would offer a similar phenomenology, such that gravity 
would be modified below a characteristic acceleration scale of the order of the one 
required in MOND. For instance if we consider the action 

S = Id 4 x^g ig ^ Gn {R ~ R 2 (Log[/])"} , (36) 

with the same assumptions on / that we did before we get that now the critical accel¬ 
eration «o also has a “running” with the mass as 


do 


G n M 



li 



48ag 

QoG n M _ 


(n—1)/2 


(37) 


To introduce some kind of scale dependence of do could be interesting since MOND 
typically gives an overestimation of the amount of visible matter at cluster scales. 

Since the expansions we have used break down for some intermediate range of 
energies/distances one should still show that the dynamics in this “non-perturbative” 
regime are consistent for some choice of / and that one recovers an acceptable matching 
between the high energy/short distance and low energy/large distance regimes. And 
then one should compute the predictions of these theories in many different situations 
for which there are experimental data to compare with before any of these models 
could be considered a viable alternative to a ACDM cosmology. This is a non-trivial 
task but it is worth undertaking it because we have seen that there are reasons to 
believe that one might explain many aspects of the cosmological and astrophysical 
observations without introducing Dark Matter in this class of theories. And, as we have 
seen, these theories also offer the unique possibility of being tested not only through 
astrophysical observations, but also through well-controlled laboratory experiments 
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where the outcome of such experiments is correlated with parameters that can be 
determined by means of cosmological and astrophysical measurements. 
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